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QUANTITATIVE ESTIMATES ON JACOBIANS 
FOR HYBRID INVERSE PROBLEMS 

GIOVANNI ALESSANDRINI AND VINCENZO NESI 


Abstract. We consider <r-harmonic mappings, that is mappings U whose components Ui 
solve a divergence structure elliptic equation div(uVui) = 0, for i = 1,..., n. We investigate 
whether, with suitably prescribed Dirichlet data, the Jacobian determinant can be bounded 
away from zero. Results of this sort are required in the treatment of the so-called hybrid 
inverse problems, and also in the field of homogenization studying bounds for the effective 
properties of composite materials. 
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1. Introduction 

The appearance of coupled physics methods has provoked a sharp change of perspective in 
inverse boundary problems. The simultaneous use of different physical modalities to inter¬ 
rogate, through exterior measurements, a body whose interior parameters are unknown has 
enabled to single out interior functionals which carry useful, and possibly stable, information 
on the parameters of interest. Such methods are also known under the name of “hybrid in¬ 
verse problems”. Notable examples are the coupling of Magnetic Resonance with Electrical 
Impedance Tomography [35], Ultrasound and Electrical Impedance Tomography m, Mag¬ 
netic Resonance and Elastography [26]. To fix ideas, let us focus on Ultrasound Modulated 
Electrical Impedance Tomography. In EIT the goal is to determine the, possibly anisotropic, 
electrical conductivity a = {or,} of a body II by repeated boundary measurements of voltage 
u\qq and current distribution a\7u ■ v with u solving the elliptic PDE 

(1.1) div(cVu) = 0, in 17. 

As is well known [35], the stability is very weak and, in fact, in the anisotropic case, also non¬ 
uniqueness occurs [30]. By combining electrical measurements with ultrasound measurements 
it is possibile to focus on a tiny spot near any point x € 17 and it has been shown by Ammari 
et al. DU that one can detect the localized energy 

(1.2) H(x) = aVu ■ V u{x). 

If one repeats the experiments with different boundary voltages, it is possible to extract the 
functionals 

(1.3) Hij{x) = (rVuj • Vuj(x). 

where u\,...u n is an array of different solutions to (jl.ip . In Monard and Bal [331 39], it 
is shown how, from such functionals, one may obtain the conductivity a in a satisfactory 
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stable fashion. The crucial point, however, is to be able to set up an array of boundary data 
(j>i ... ,4> n and corresponding solutions u\, ... u n in such a way that the functionals Hij are 
non degenerate. 

In other words, calling U : fl —> M n , the mapping U = (u\, ..., u n ), which we shall designate 
“cr-harmonic” mapping, it is required that the Jacobian determinant 

det DU 

does not vanish. And, furthermore, for the purpose of stability, a quantitative lower bound 
would be needed. 

This is the main question that we wish to address in this note, which essentially stays behind 
all coupled physics problems mentioned above, and other inverse problems as well. The same 
issue showed up, for instance, in the field of groundwater transmissivity detection [43]. 

This kind of questions also arises in the branch of the homogenization theory which studies 
effective properties of composite materials. We give a brief outline here. 

Indeed, the positivity of Jacobians of injective cr-harmonic mappings has attracted attention 
in several applications. In two dimensions, the first application of this positivity has been 
given in [42]. The long standing problem of improving the so-called Hashin-Shtrikman bounds 
[25] for the effective conductivity of composite materials was addressed in that paper. The 
method used was based on ideas of Murat and Tartar m and Tartar [46], a reference not 
easy to find. We refer to m for a more complete treatment. 

The bottom line is as follows. The question of interest, in the simplest not yet known at 
that time, case is the following. Three numbers 0 < o\ < (72 < CT 3 , representing the con¬ 
ductivity of three isotropic materials, called the phases, three “volume fractions” pi,P 2 and 

P 3 , summing up to 1 and representing the area proportions of the phases and a 2 x 2 matrix 

3 

A, parametrizing the affine boundary data, are given. Assume that cr = E Xi( x ) a i where 

2=1 

Xi represents the characteristic function of the set where o is equal to cq times the identity 
matrix and py f n Xi( x )dx = Pi,i = 1,2,3. Then one aims to determine a bound from below 
for the following quantity 

(1.4) F(A) = inf inf —f Trace[(Ut/o(x) + A) t ct{x){DUq{x) + A)\dx . 
xi,X2,X3 UoeW ^ 2 ^) |S2| J n 

The overall problem is non linear and actually is linked with the notion of quasi-convexity. 
More precisely, it computes the quasiconvexification at the matrix A, of a non-convex function 
of DU, with U — Ax € Wq’ 2 (Q). This function turns out to be the minimum of three quadratic 
functions, as shown by Kohn and Strang in [29]. However, for our purposes, it is important 
to note that the infimum over Uq in (11.41) is attained exactly when U(x ) = Uq(x) + Ax is 
the cr-harmonic mapping with affine boundary data given by U = Ax on dfl. Optimal lower 
bound for (11.41) were found by Kohn and Strang exactly exploiting the connection with the 
optimal bound for effective conductivity found by Murat and Tartar and, later, by Cherkaev 
and Lurie [33]. The optimality is restricted to the case when only two isotropic phases are 
present that is, only two materials are “mixed”. For three or more phases, the methods 
based on compensated compactness gave suboptimal bounds. In this specific context, the 
compensated compactness method uses simply the constraint that the Jacobian determinant 
of the matrix DU is a null-lagrangian. 
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The classic strategy gives the so-called “Wiener bound”, that is the harmonic mean bound. 
It is obtained considering the test fields B in the class 


ix = A > . 


One obtains 


(1.5) 


So := € L 2 (Q) : -L f n B (x)d. 

F(A) > F 0 (A ) := inf f Trac e[B(x) r a{x)B(x)]d. 

B&Bo |XZ| j n 


Trace 


A 1 


N 


f a 1 (x)d: 
Jn 


-1 
x ) A 


Tartar’s ideas, based on compensated compactness, in this simplified context lead to an 
improved bound (called the “translation bound” by G. W. Milton) obtained by considering 
the new test field in the class 

B\ := 1 B G L 2 (G) : —— f B(x)dx = A, —— f det B(x)dx = det A 

[ l“| Jn l“l Jn 

One obtains 



F(A) > Fi(A) := inf -i- f Trac e[B(x) T a(x)B(x)\dx. 

The computation is more involved than (11.51) . This technique, however, gives an optimal 
answer in two dimensions, when specialized to the case of two-phase isotropic materials. 
When one deals with more than two phases this approach is no longer optimal. The results 
in m have the following corollary. Set 


£>2 := | B € L 2 (G) : —— f B(x)dx = A,—— f det B(x)dx = det A det B > 0 , a.e. in O 

l |»‘l Jn l“l Jn 


One has 


(1.6) F(A) > F 2 (A) := inf -j— f Trac e[B(x) T (r(x)B(x)]dx. 

bgb 2 |i2| Jq 


In fact, in [32], it is proved that F 2 (A) > F\ (T), as soon as one deals with more that two 
isotropic phases, for suitable choices of the given parameters pi and matrices A. Later, new 
optimal microgeometries were found for multiphase materials in [2] and, using again the 
positivity of the Jacobian determinant, it was possible to prove their optimality according 
to a stricter criterion, see [3]- The key is exactly the universal bound given on the Jacobian 
determinant, which, in this context reads as the inequality det A det B > 0 in (11.61) . In this 
context, it is highly desirable not to have any constraint on the regularity of the interfaces 
between phases. When a is non-symmetric, applications to composites have been given, 
for instance, in the context of the classic Hall effect by Briane and Milton [15l 16] . Other 
applications have considered the problem of determining which electric fields are realizable 
by Briane, Milton, and Treibergs[18]. On the other hand one would like to have similar 
improvements in higher dimensions. Briane and Nesi [T9] studied the case of laminates of 
high rank showing that, for these special microgeometries the positivity of the determinant 
of the “Jacobians” of the corrector matrix holds in any dimension. To explain the result 
in detail would require too long a digression. However, roughly speaking, one could expect 
that in higher dimensions, even for discontinuous a one could hope for the positivity of the 
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Jacobian determinant if one makes assumptions on the “microgeometry”. On the other hand, 
even in the very restricted setting of periodic boundary conditions, particularly adapted to 
composites, and even under the assumption of dealing with only two isotropic phases, there 
is no hope to control the sign of the Jacobian determinant of cr-harmonic mappings without 
further assumptions on the nature of the interfaces. One explicit example was provided by 
Briane, Milton and Nesi m- 

We now go back to the precise subject of the present paper. We pose the following problem. 

Problem 1. 

Can we find Dirichlet data 

(1.7) $ = (&,- ■■ An) :dn^R n 

such that the corresponding solution mapping U = ,u n ) is such that det DU is 

bounded away from zero independently of the conductivity cr? 

Note that, in this context, it is essential that the choice of the boundary data is independent of 
cr, because a is the real unknown of the original inverse problem. As is easily understandable, 
some a-priori assumptions on cr, such as ellipticity, and some kind of regularity shall be needed. 

Problem 1 has a different phenomenology depending on the space dimension. When n = 2 
the issue is more or less completely understood, whereas when n = 3 or higher, various kinds 
of pathologies show up. A review of such pathologies and a discussion of the open issues when 
n > 3 shall be the object of Section 01 

The principal aim of this note is to provide, when n = 2, a quantitative lower bound on the 
Jacobian determinant under essentially minimal regularity assumptions. This is the content 
of our main Theorem 12.111 which is the new contribution of this paper to this subject. 

We start reviewing the main known results in dimension n = 2. It was proved in Bauman et 
al m that, if cr is Holder continuous, H has C 1,a boundary and is a C 1,a diffeomorphism 
onto the boundary of a convex domain, then det DU > 0 everywhere. Note that in m, 
only symmetric matrices cr were explicitly considered, however, in view of classical results 
on two dimensional elliptic first order systems with Holder coefficients see, for instance, m 
Appendix and also m Proposition 5.1, the result extends as well to the non-symmetric 
case. On the other hand, the present authors [7], proved that when a is merely L°° and T is a 
homeomorphism onto the boundary of a convex domain, then det HU > 0 almost everywhere. 
In fact it was proved that, for every locally invertible, sense preserving, cr-harmonic mapping 
U one has 

(1.8) log det DUeBMO 
and, subsequently [8], this result was improved to 

(1.9) det DU € Aoo 
that is the class of Muckenhoupt weights [21]. 

We recall that for purely harmonic mappings, Lewy’s Theorem [33], states that for two- 
dimensional harmonic homeomorphisms, the Jacobian determinant cannot vanish at interior 
points. Hence, when n = 2, harmonic homeomorphisms are, indeed, diffeomorphisms. How¬ 
ever the Jacobian determinant may vanish at boundary points. 
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It is also worth mentioning that the convexity assumption on the target of the boundary 
mapping <b is sharp, Choquet [ 2 Q], Alessandrini and Nesi [9], if one wishes to have a condition 
expressed merely on the “shape” of the target and not on its parametrization. 

Conversely, note that when no regularity is assumed on a, the essential infimum of det DU 
on compact subsets of D might indeed be zero. In Section [3] an example, based on a well- 
known one by Meyers, is illustrated. 

In the next Section [2] we shall prove a quantitative version of the result in m- The starting 
point relies on prescribing some quantitative assumption on the boundary data <h, when 
viewed as a parametrization of the boundary of the convex target, see Definitions 12.11 12.21 
12.61 The subsequent step consists on a quantitative lower bound of the modulus of the gradient 
of a scalar solution to equation (11.11) , Theorem 12.91 This estimate may be interesting on its 
own. Finally we state and prove our main result, Theorem 12.111 

2. The quantitative bounds. 

Let f : M — > M be a T-periodic C 1 function. Let u : [0, oo) — > [0, oo) be a continuous strictly 
increasing function such that cj( 0 ) = 0 . 

Definition 2.1. Given m. M € M, m < M, we say that f is quantitatively unimodal if there 
exists numbers t± < t- 2 < £3 < £4 < f 1 + T such that 



fif) 

= m 

t € [ti,t 2 ] , 

( 2 . 1 ) 

fit) 

fit] 

-fit) 

= M 

t € [t 3 ,U] , 

> min{w(t — t 2 ), w(t 3 — t)}, 

t <E [* 2 ,^ 3 ] , 


> min{u;(t — + T — f)}, 

t £ [£ 4 , t\ + T] 


In the sequel we will refer to the quadruple {T,m, M,u} as to the “character of unimodality” 

off. 

The concept of unimodality, but not this terminology, first appears in Kneser [28] , when 
he proved Rado’s conjecture [44] concerning the case of “purely” harmonic mappings. The 
terminology “unimodality” was introduced in this context by Leonetti and Nesi [32], following 
the work of Alessandrini and Magnanini [ 6 ]. A different terminology (almost two-to-one 
functions) has also been used for the same concept, Nachman, Tamasan and Timonov[41J. 
Let T C l 2 be a simple closed curve parametrized by a T-periodic C 1 mapping 

(2.2) $ : R -A M 2 

in such a way that ^| [0T) is one-to-one. 

Definition 2.2. We say that T is quantitatively convex if for every £ e R 2 , |£| = 1 the 
function 

ft = $ • £ 

is quantitatively unimodal and its character of unimodality is given by {T, mg,Mg,u;} with 
mg, Mg such that Mg — mg > D, for a given D > 0. 

In the sequel we will refer to the triple {T,D,tv} as to the “character of convexity” ofT. 

Remark 2.3. If T is quantitatively convex then it is convex, that is, it is the boundary of 
a convex set G. In fact each tangent line to T turns out to be a support line for G. The 
following Lemma provides a sufficient condition for quantitative convexity. Roughly speaking, 
it says that if T is an appropriately parametrized C 2 simple closed curve with strictly positive 
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curvature, then it is quantitatively convex in the sense of Definition \2.2l and the character of 
convexity can he computed in terms of the parametrization. Here, for the sake of simplicity, 
we have chosen the arc-length parametrization, because the main purpose of this Lemma is 
to provide a variety of examples, but we emphasize that in general, the character of convexity 
does depend on the parametrization of the curve and not only on its image. 


We convene to denote by J the matrix representing the counterclockwise rotation of 90 degrees 



Lemma 2.4. Let T be such that $£C 2 and: 


(2.3) 


i) |<f>'| = 1 , 

ii) 0 < k < 4*" • < K , 


then T is quantitatively convex with character {|F|, ^ t}. 


Proof. Condition i ) of Lemma 12.41 implies that &(t) = e lsd \ 0 < t < T and we may assume 
that s(0) = 0. Without loss of generality we assume that 4> is orientation preserving. Then, 
by condition ii) of Lemma 12.41 one has 0 < n < s'(t) < K. Picking, w.l.o.g., £ = e^, 

M t )~ ( l ) d 0 ) = f sin (s(r))dr. 

Jo 

The function s(t) ranges over the whole interval [0, 2tt\, picking t n such that s(t 7r ) = 7 r, we 
have 

M t ~ = ft ds (r ) > £ fl sin (s)ds = ± , 


= sin(s(i)), 

s(t) = f^s'(r)dT > nt, 

M*) — ^rf t > 0 ^ ^ f ■ 

Thus we may pick D = -^ and oj(t) = ^ t , t > 0. □ 

We shall consider fl a bounded simply connected domain in M 2 with C 1,a boundary. In 
order to make precise the quantitative character of such regularity we introduce the following 
definition. 

Definition 2.5. A domain 14 C K 2 is said to be of class C 1,a with constants pq,Mq, positive 
and Holder exponent a € (0,1], if for any P E dLl, there exist a rigid change of coordinates 
such that P = 0 and we have 


(2.4) 

where if : 

(2.5) 
and also 

( 2 . 6 ) 


[ POi Po] 


n n B po ( 0) ={i£ B po (0) : x 2 > if (aq)}, 
is a C 1,a function satisfying 

if{ 0) = ^'(o) = o 


1 +a 


Ia°°([-Po,Po]) + P 0 11 if' 11 L°° ([— po ,Po] ) + Po ' SU P 

x,x'e[-p 0 ,p 0 \ 

x^x' 


\lf\x) -lf\x')\ 


\x — X 


/let 


< MqPq . 
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Definition 2.6. Given a C 1,a (dFl\ M) function (j), we shall say that it is quantitatively uni- 
modal, if considering the arclength parametrization of dLl, x = x(s), 0 < s < T = |<9D|, the 
periodic extension of the function [0, T] 3s-> (f>(x(s)) is quantitatively unimodal with char¬ 
acter {T,m, M,u}. For such a function <f>, we introduce the following closed arcs, possibly 
collapsing to a single point: 

, . r min = {i£ dLl : (j) = m} , 

[ } r max = {x £ dCl : (j) = M} . 

Accordingly, a mapping <1? G C 1,a (dLl; M 2 ) shall be said quantitatively convex with character 
{T,D,iw} if the periodic extension of$>(x(s)) fulfils the conditions of Dehnition \2.'A . 


Let us consider a = 1,2 a, not necessarily symmetric, matrix of coefficients : Ll —>• K 

satisfying the ellipticity condition 

> K - 1 |£| 2 , for every £ G M 2 , 

’ > K ~ l |£| 2 , for every ^ G M 2 

for given positive constant K, and also 

(2.9) \aij(x) — aij(x')\ < E\x — x'\ a , Vx, x' Gfl, 


for given a, 0 < a < 1 and E > 0. 

We shall consider the W 1 , 2 (Q) solution u to the Dirichlet problem 


( 2 . 10 ) 


( div(cVu) = 0 in , 

1 u = (f> on dFt. 


We recall that, in view of the classical regularity theory, u in fact belongs to C' 1 , ^(fl), for 
some (3 < a and its norm is dominated by the C' 1,a -norm of <f>, modulo a constant which only 
depends on pq,Mq,K and E, with po, Mq as in Definition 12.51 


Lemma 2.7. Let 4 > : dLl -A M be quantitatively unimodal with character {\dLl\,m, M,uj} and 
assume that 

( 2 . 11 ) \-^(j)(x(s)) --^(j}(x(s'))\ < E\s - s'\ a , Vs, s' € [0, |<9f2|]. 

Then there exist k, 6 only depending on the character of unimodality (see Deftnitions \2. 1[ \2.b}) 
and on a, E, such that if 

(2.12) x G and dist(x, r m ; n U r max ) < 5 , 
then 


(2.13) 


|Vu(x)| > K. 


Proof. Up to a C 1,a diffeomorphism, with constants only depending on po, Mq and |5D|, we 
may assume that D = i?i( 0 ). 

It is well known that in such new coordinates u solves a new Dirichlet problem of type (12.1011 
with a new matrix of coefficients and new boundary data that, however, satisfy analogous 
assumptions with constants and parameters only depending on the same a-priori data. For 
the sake of not to overburn the notation we stick to the one of (| 2 . 10 l) . 

By the C 1,/3 regularity of u, if dist(x, r max ) < 77 , then u(x) > M — Crj with C > 0 only 
depending on the a-priori data. 

Let us pick 77 such that 


M — Cp > 


M — m 


2 
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Hence, by Harnack’s inequality [24] , 

u(x) — m>C AI f m > 0, for every x £ B 1-^(0). 

Here C only depends on the a-priori data. By the version of the Hopf Lemma due to Finn and 
Gilbarg [23] Lemma 7, which applies to equations in divergence form, and Holder continuous 
a, we obtain 

(2.14) |Vit(x)| > no > 0 , Vx € r min , 

with Ko only depending on the a-priori data. 

By C 1 ' 13 regularity we have 

|Vrt(x)| > kq — C5 13 Vx € fi such that dist(x,r m ; n ) < 5. 

Picking S such that C5 3 < we obtain 


(2.15) |Vrt(x)| > > 0 , if dist(x, r m ; n ) < 6 . 

A symmetrical result applies in the neighborhood of r max . □ 

Lemma 2.8. Under the same assumptions as in Lemma \2. 1\ there exists r > 0 such that 

(2.16) |Vu(x)| > L > 0 , Vx G O, dist(x, dQ) < r . 

Here L and r are positive and only depend on the a-priori data. 

Proof. If we pick x £ dQ, and write x = x(s) such that dist(x, r m i n U r max ) > 6, we have 


|Vu(x(s) • a/(s)| 




> cv{6). 


By C 1,/3 regularity 

|Vu(x)| > min{K, w(4)} — Cr 13 , Vx G such that dist(x, <9f2) < r . 
Picking r such that 

Cr 3 = i min{K, w(4)} , 

the thesis follows. 


□ 


Theorem 2.9. Let Ll be a simply connected domain, C 1,a -regular with constants {pq,Mq\ 
(see Definition IH.,51) . Let (f) : dil —>■ R be quantitatively unimodal with given character 
{\dLl\,m,M,uj}(see Definitions ] 2 . il [OP and let it satisfy the Holder condition (12.111) . Let 
a = {cr ? j (x)},; J= i i 2 satisfy the ellipticity condition (12.81) and the Holder bound (12.91) . Let 
u £ be the solution of the Dirichlet problem (12.101) . 

Then there exists C > 0, only depending on the a-priori data as above, such that 

(2.17) |Vrt(x)| > C > 0 , for every x £ D . 

Remark 2.10. Under stronger regularity assumptions, in particular assuming that a is Lip- 
schitz continuous, a similar result was proven already in [4] ; Theorem 3.2. 

Proof. As is well-known, there exists u £ called the stream function associated to 

u, which satisfies 

(2.18) Vu = J(jVu everywhere in II, J = 
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Using complex notation z = x i + ix 2 , f = u + iu, the system (12.181) can be rewritten as 
(2.19) fz = Vfz + vfz inH, 


where, the so called complex dilatations /i, u are given by 


( 2 . 20 ) 


_ 0-22—oil— ^(o'l2+g'2l) _ 1 —det cr+i(o\2— an ) 

1 1 +Tr cr+det a ’ 1 +Tr er+det a 


and satisfy the following ellipticity condition 


( 2 . 21 ) 


H + M < 


K - 1 
K + 1 ’ 


and, being cr Holder continuous, also ji and v satisfy an analogous Holder bound. 

In P2], it is proven that / is a C 1 ’^ diffeomorphism of onto /(H). The lower bound obtained 
in Lemma 12.81 implies that, setting 


= {x G fl : dist(x, dfl) > r}, 

/ : fl\H r —>■ C, is a bilipschitz homeomorphism with constants only depending on the a-priori 
data. We have identified C with R 2 in the canonical way. 

Hence /(H) is also a domain with constants controlled by the a-priori data. Note that 
also |<9(/(H))[ is controlled. 

Let us denote g = w £ C. A straightforward calculation gives 

( 2 . 22 ) = - v{g)g w - g{g)g^ ■ 


In other words g satisfies a Beltrami equation whose coefficients satisfy uniform ellipticity 
and Holder continuity, with constants only depending on the a-priori data. 

By standard interior regularity estimates, g w is bounded in /(H r ). Using (12.221) . we have 

(2.23) \g w \ 2 ~ \gui\ 2 < C 2 in /(H r ), 


which can be rewritten as 

\fw\ 2 - \fw\ 2 > C~ 2 in n r , 

which in turn implies 

(2.24) | Vu[ > C- 1 in H r . 

Hence, in combination with Lemma 12.81 the thesis follows. □ 


Theorem 2.11. Let H and a be as in Theorem 12. hi Let <h = (/q,^) : dLl — > R 2 be 
quantitatively convex, see Definitions \2 . HI POl with character {\dLl\, D,uj}. LetU = ( 111 , 112 ) € 
W 1 ' 2 (H;R 2 ) solve 

/ 2 25 n f di v(aVui) = 0 in H , 

^ ' ) Ui = 4>i on dVt. 

There exists C > 0 only depending on the a-priori data such that 

U-.n^ U(Q) C R 2 

is a C 1 ^ diffeomorphism and 

(2.26) det DU > C 2 > 0 in H. 
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Remark 2.12. A similar result, under slightly more restrictive hypotheses, has been recently 
proved by G. S. Alberti |Tj. In fact the approach in [lj is based on estimates in |4j and [5] 
which require the Lipschitz regularity of a. Conversely, under somewhat different regularity 
assumptions, quantitative upper bounds on the so-called dilatation of a a-harmonic mapping 
U, that is the quotient 

Tmce(DU T DU ) 

2 det DU ’ 

have been recently studied in [TU] , Theorems 3.1, 3-4 ■ 

Proof. In m, it is shown that U is a orientation preserving diffeomorphism. The lower bound 
(12.261) remains to be proven. For any £ G K 2 , |£| 2 = 1, we may apply Theorem 12.91 to 

U£ = U-Z 

and obtain 

(2.27) \(DU)f\ = |Vu 5 | > C > 0 in H. 

Or equivalently 

(2.28) \DU T DU£-t\ = |Vu s | 2 > C 2 > 0 in U, 

that is the eigenvalues Ai(x) and \- 2 (x) of the symmetric matrix DU T (x)DU(x ) are uniformly 
bounded from below: 

\i(x) > C 2 > 0, i = 1,2, Vxef2. 

Therefore 

(2.29) (det DU) 2 = \i(x)\ 2 (x) > C 4 > 0 , everywhere in 0. 

Since U is sense preserving, one has det D U > C 2 > 0 everywhere in 0. □ 

Remark 2.13. Theorem, \2. 1 11 has some feature in common with the results in \ 9]. In the 
latter paper the authors consider harmonic mappings which are extensions of given Dirichlet 
data aiming for univalent solutions. A characterization is given, for the case when U is a 
diffeomorphism up to the boundary of Q, in terms of the value of the Jacobian determinant 
on the boundary, so implicitly imposing constraints on the parametrization of the boundary 
of the image. One may wonder whether an assumption just on the shape of the target may 
suffice. This is not the case even in the purely harmonic case. Indeed one may exhibit a 
sequence U n of sense preserving, injective, harmonic mappings of the unit disk onto itself, 
fixing U n (0) = 0, such that detDU n (0 ) —> 0 as n —>• +oo. The convergence holds uniformly 
on compact subsets of the unit disk. The limit harmonic mapping in not univalent. See [22j . 
Section 4-1- 


3. Discontinuous coefficients. An example. 


We elaborate on a well-known example by Meyers [37|. See also Leonetti and Nesi [32 ] for 
an application in a related context. For a fixed a > 0 we consider the symmetric matrix of 
coefficients 


(j(x) = 


( a 1 x\+ax\ (o 1 — a)xiX 2 \ 

-!T^-n- -ri-n- 



(3.1) 
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Is is a straightforward matter to check that its entries belong to L°° and that a has eigenvalues 
a and a -1 . Therefore a satisfies the uniform ellipticity condition (|2.8j) with ellipticity constant 

K = max{a, a -1 } , 

and a is discontinuous at (0,0) (and only at (0,0), when a / 1). Let us denote 

ui(x) = , 

U2 ( x ) = \ x \°‘~ 1 X2 ■ 


A direct calculation shows that tq £ W 1,2 (Bi(0)), i = 1,2 and that they solve the Dirichlet 
problem 

div(fjVrtj) = 0 in L?i(0), 

Ui = Xi on dB i(0). 

Note also that f = u\ +iu 2 is a quasiconformal mapping of -Bi(O) onto itself and it solves the 
Beltrami equation 

h=^rifz- 

a + 1 z 

Setting U = (iq , 112 ), we compute 

detDU = \f z \ 2 -\tt 2 = a\z\ 2 ( a ~ 1 K 


Therefore det DU vanishes at (0,0) when a > 1, whereas, when a £ (0,1), it diverges as 
z -> 0. 


4. Mappings in higher dimensions. Examples and open problems. 

The interior lower bound on det DU obtained in Theorem 12.111 has been achieved by 
methods which are intrinsically two-dimensional (the Beltrami equation). Only part of the 
result can be extended to higher dimensions. 

For instance, with minor adaptations of the method developed in the Section [2j one can argue 
as follows. 

Consider 17 C M n , a bounded domain diffeomorphic to a ball of class C l,a and with constants 
Po,Mq defined with the obvious slight adaptations of Definition 12.51 

Let a = {cr t j} l j = i 2 be the matrix of coefficients and let it satisfy uniform ellipticity with 
constant K as in (12.811 and Holder continuity like in (12.9|) . 

Let G C M n be a convex body whose boundary T is C 2 and having at each point principal 
curvatures bounded from below by k > 0. 

Let $ = (cj)i,(f> 2 , ■ ■ ■ , (f>n) '■ <912 -A T be an orientation preserving diffeomorphism such that 
<4>, <L» —1 are C 1,a with constant E. Let U = (u\,U 2 ,--- , u n ) £ fT 1,2 (12; M n ) be the weak 
solution to 

{ div(crV'Uj) = 0 in 12, 

Ui = 4>i on 12 , 

i = 1 , 2 • ■ • , n . 

Then, by the same arguments used in Section [2l we obtain. 

Theorem 4.1. Under the above stated assumptions, there exists p > 0 and Q > 0 such that 
U is a diffeomorphism ofQ\Q p onto a neighborhood ofT, within G and we have 

det DU > Q in 12\12p . 
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We omit the proof. 

When n > 3, there is no chance, under the kind of hypotheses stated above, to obtain a 
global lower bound on det DU. Evidence comes from a sequence of counterexamples that 
have been produced in a wide time span. A first illuminating example goes back to Wood [48] 
and has the amazing feature of being totally explicit. Wood displayed the following harmonic 
polynomial mapping from M 3 onto M 3 : 

U(x 1 ,x 2 ,x 3 ) = (ui,u 2 ,u 3 ) = (x? - 3xi x\ + x 2 x 3 , x 2 - 3xix 3 ,x 3 ) 

that is U is a homeomorphism, but not a diffeomorphism because det DU = 0 on the plane 

{xi = 0 }. 

Later Melas [36] provided an example of a three dimensional harmonic homeomorphism 
U : B i(0) —> B i(0) such that det DU(0) = 0. Subsequently, Laugesen [31], showed that there 
exists homeomorphisms <f> : dB i(0) —>■ dB i(0) which are arbitrarily close to the identity in 
the sup-norm such that the mapping U = (u\, u 2 , • • • , u n ) solving 


Aui = 0 in B\ (0), 

U = $ on 0£i(O), 

i = 1, 2 • • • , n . 


is such that det DU changes its sign somewhere inside -Bi(O). Such examples are especially 
striking because they show that in dimension n > 3 it seems difficult to find a universal rule to 
select Dirichlet data in such a way that the corresponding harmonic (or cr-harmonic mapping 
is invertible at a topological level (because it may reverse orientation!) and not only as a 
differentiable mapping. 

One may wonder whether changing the topology of the boundary data may help. In the 
periodic case, obviously the harmonic functions are linear and this might have left the hope 
that, for variable coefficients, the periodic case may be better that the generic Dirichlet 
problem. However this is not the case. In m it was proved that, in dimension three, 
one can find a matrix cr taking only two values, proportional to the identity matrix, and a 
periodic arrangement with a smooth interface, but such that the corresponding solution U of 
the cell problem also reverses the orientation. The Jacobian determinant changes its sign in 
the interior of the (unit) cube of periodicity. 

If, from the above examples, it seems that few chances are left of finding a universal criterion 
by which choosing Dirichlet data such that, for each a (although smooth) the corresponding 
cr-harmonic mapping U has nondegenerate Jacobian, then a more reasonable goal would be 
to find a way to control, in term of the Dirichlet data, the set of points where the Jacobian 
may degenerate and possibly evaluate the vanishing rate at such points of degeneration. 

This appears as a completely open problem, not at all easy as the following example by 
Jin and Kazdan [27] shows. Let a € C'°°(M;M) and set 


(4.1) 


with 



( 

1 

a(x 3 ) 

0 


a(x) = 

a(x 3 ) 

1 

0 


V 

0 

0 

b( x 3) 

a(x 3 ) 

= 0 

for 

x 3 < 0, 


a(x 3 ) 

£ (0, ao) 

for 

x 3 > 0 

with 

b(x 3 ) 

l 

1— d z {xz) 

for 

x 3 €R. 



°0 € ( 0 , 1 ), 
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We set 

(4.2) U ( X ) = (X!,X 2 , -X!X 2 + <j>{x 3 )) , 

where f is chosen in such a way that 

(i bqyy — 2a = 0 , .T3 G I, 


4>(x 3 )=0, x 3 < 0 . 

It turns out that $ > 0 for x 3 > 0 and consequently 

£ 3<S: 

This means that the Jacobian determinant of a er-harmonic mapping does not fulhll the 
property of unique continuation (whereas this is the case for \DU\ 2 = Trac e(DU T DU)). 
Hence the evaluation of the zero set of det DU from boundary data might be troublesome. 

Remark 4.2. The above example has some striking features. First of all note also that, 
letting ao \ 0, we can make a as close as we want to the identity matrix. Moreover U 
converges, uniformly on each compact subset of M 3 , to the harmonic polynomial mapping 
Uq{x) = {xi,x 2 , -xix 2 ). 

We conclude by noticing that a limiting case of the above construction yields an example 
with a discontinuous, two-phase, u which is remarkable as well. 

As before we pose 

f 1 a(x 3 ) 0 

(4.3) a(x) = I a(x 3 ) 1 0 

\ 0 0 b(x 3 ) 

where now 


0 ( 2 : 3 ) = 0 
0 ( 2 : 3 ) = a 0 

K*3) = i-J^) for x 3 €R. 
That is a is piecewise constant, namely 


for x 3 < 0 , 

for x 3 > 0 with ao € (0,1), 


( 1 

0 

0 \ 

( 1 

ao 

° \ 


0 

1 

0 

when x 3 < 0 and a(x) = a o 

1 

0 

when x 3 > 0 

\o 

0 

1 / 

{ 0 

0 

l— a o / 



Again we pose 

(4.4) U(x) = (xi,x 2 , -xix 2 + <t>{x 3 )) , 

where now f is given by 

4>(x 3 ) = o 0 (l - al)xl , 2:3 > 0 , 


f(x 3 ) = 0 , 2:3 < 0 . 

We obtain that U is a a-harmonic mapping with C 1,1 regularity and, analogously to the 
previous example, it satifies 


det DU = 


J 2ao(l — Oq) x 3 > 0 , for x 3 > 0 , 


0 , 


for x 3 < 0 . 
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